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CHAPTER FIVE 

MULTIPLE SCATTERING BY AN EXTENDED REGION OF 

INHOMOGENEITIES 

In this chapter we shall treat two types of extended inhomogeneities: (i) periodic and 
(ii) random. 

1 Waves in a periodic medium 
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Propagation of light or sound wave is of long standing interest in several branches of 
basic and appied physics, from old disciplines such as x ray diffraction in crytallography, 
to the modern science of photonic crystals. Many problems in natural environment also 
involve wave propagaion in periodic media. For example, nearly periodic sand bars are 
freqeuently found in shallow seas outside the surf zone; their presence changes the wave 
climate near the coast. The technology of remote-sensing, either by underwater sound 
or by radio waves from a satellite, depends on our understanding of scattering by the 
wavy sea surface. 
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Figure 1: Bragg resonance due to contructive interference 

In a periodic medium, the pheonomenon of Bragg resonance holds a special position. 
Let us explain it for the one-dimensional case of monochromatic water waves passing over 
a long stretch of parallel bars on an otherwise horizontal seabed. In nature, the bars are 
usually of much smaller amplitude than the wave depth. Most waves can pass over them 
without singificant reflection, except when the wavelength is an integral multiple of the 
bar period. In Figurel, the special case where A = 2\bar is sketched. Upon encountering 
each bar crest, every incident wave crest is mostly transmitted and slightly reflected. 
At a given bar crest, the height of the reflected wave crest is the sum of infintely many 
left-going crests. We called it the n— th crest if its reflection is originated at the n-th 
bar crest on the right. Because of the 1 : 2 ratio, each crest has traveled the distance 
of 2n\bar = nX since its first passage over the bar crest in question. As a consequence, 
they are all in phase with one another, hence the sum of the reflected wave intensity is 
very high. 

We treat first waves in an elastic laminate of infinite extent. Borrowing the existing 
knowledge in solid-state physics, we discuss the relation between Bragg resonance and 
band gaps in the dispersion relation. We then treat the more practical case of finite 
extent in order to study the tranmission and reflection. The asymptotic method of 
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multiple scales is applied to derive the coupled-mode equations for the slowly varying 
envelopes of incident and scattered waves. 

The physics discussed here is representative of two and three dimensional period 
media, which would call for a much more elaborate theoretical analysis. 



2 Waves in an infinitely thick laminate 

We begin with 

for the longitudinal displacment u{x,t). We shall assume that the elasticity is peri- 
odically modulated about the mean value Eq. Defining E{x) — Eq£ where £^ is a 
dimensionless periodic function of x. Equation (2.1) can be rewritten as 

where C — -^/Eq/ p. 

Consider a monochromatic wave in an elastic laminate so that 

u{x,t)^lk{U{x)e-''^^) (2.3) 

then 

IK} 4^-° P-^) 

What sort of wave can exist with period L — Na (Born- von Karman condition)? Math- 
ematically this is a homogeneous BVP, hence only eigensolutions exist for certain value 
of the eigen frequency uj. In the limiting case of uniform medium £ = \, the solution is 
clearly 

U - Ve'^'' (2.5) 

where uj — Ck. 

To solve the eigenvalue problem, Bloch (and Floquet) found the following theorem. 
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2.1 1-D Bloch's theorem 



Bloch' theorem: For the one-dimensional problem governed by the ODE 

4"- ( £{x)^^] + a^U = 0, where £{x) = £{x + a) and a = (2.6) 
ax \ ax J C 

subject to the condition that U{x) is periodic over the period L = Na, the solution is 
of the form: 

U{x) = U{x)e'^'' (2.7) 

where U{x) — U{x + a) is periodic in x and k is an eigenvalue. 

We follow the reasoning in Ashcroft & Mermin in their more general proof for 3D 
quantum mechanics. 

Since E{x) is periodic with period a we can Fourier expand 

£{x) = ^FKe'^'' (2.8) 

K 

where K is the short-hand notation for — 2n7r/a, with n = 0, ±1, ±2, ±3, ... ± oo 

and Fk is a function of K. Note that 

^iKa ^ ^2nn ^ ^ ^2.9) 

Let the general periodic soluton be 

U{x) ^J^^i^"^"" (2-10) 

where the expansion coefficient is a function of q. We substitute (2.8) and (2.10) into 
(2.6), to get 

a' J2 Ke'"' - E E + K)FKV,e'^'^+''^^ = 

q q K 

Let us change the summation variable in the double series by letting q + K — q' , or, 
q'^k-K: 

a' J2 Vqe'''' - E - K)q'FKV,,^Ke"''^ = 

q q' K 

or, 

„2ta \ V„ r^\„TP M ^iqx 



E 



V, - E(<! - K)qFKY,.K 



a 

K 



Becasue of the Born-von Karman condition, it can be shown that for different q, e*^^ 
are orthogonal to each other, so that 



X - Y.^<1 - K)qFKV,^K = (2.11) 

K 

Let us change the dummy symbol K to K' in the sum, and replace 

q^k-K (2.12) 
a''Vk-K - $^(fc -K-K'){k- K)FK'Vk-K-K' = (2.13) 

K' 

Let us change again: K" — K + K' in the sum, then 

a^Vk-K - J](A; - K"){k - K)FK"-KVk-K" = (2.14) 

K" 

Because K = Kq, ±Ki, ±K2, this is an infinite set of simultaneous equations for 

For a fixed k, we have by definition (2.12), q — ...k + K2, k + Ki, k,k — Ki, k — K2, k — 
K3, .... The corresponding coefficient solutions are: ...Vk+K2, Vk+Ki,Vk, Vk-Ki, Vk-K2, ■■■■ 
Thus the wave solution is 

Uk{x) = Y,Vk-Ke'^'~''^^ (2.15) 

K 

Since for a fixed k 

Uk{x) = Vk-Ke''^^-''^' = I ^fc-i^e-*^^ I 6^^=^ (2.16) 



K \ K 



Denoting 



we get 



Clearly 



Uk{x) = YVk_Ke-'''^ (2.17) 



K 



Uk{x)=Uk{x)e''''' (2.18) 



Uk{x) = ^fe-i^e-^^- = Uk{x + a)^Y 14-Ke-'^"e-^^« (2.19) 

K K 

because of (2.9). Bloch's theorem is proven for the one dimensional problem. 
We now use Bloch's theorem to solve the eigenvalue problem. 



2.2 The dispersion relation- the eigenvalue condition 

Consider a monochromatic wave in an elastic laminate so that 

u{x,t) = ^{U{x)e-'^^) (2.20) 

then 

where — Eq/p and S is periodic in x with period L. In in general £ can be expanded 
as a Fourier series 

oo 

S{x)= e^e''"^" (2.22) 

m=— oo 

where K — ^ is the wavenumber of the fundamental harmonic^. We can assume that 
S consists of the mean and the fundamental harmonic and rewrite it as 

£{x) + eme'"^"^^ (2.23) 

Applying the Floquet-Bloch theorem we assume the solution to be of the following form: 

U{x) = 6^^=^ J2 yn{x)e'''^^ (2.24) 

n 

Substituting (2.23) and (2.24) and in (2.21), and noting that 

e±^^-^ = ^ i(^ + nX)Ke^['+("^^)^l" = ^ + (n ^ l)i^]KTie'[^+"^]^ 
we get 

5Z - + ^^)') Vne'^^^^''> - emVnik + nK){k + {n + m)ir)e*(*^+("+"^»^ 

n ^ ' n,m 

After changing the summation indices from n + m to n', we find 
J] - (A; + nK)A yj^^^'^'^'^'' - ^ e^V;,_^(/t + (n' - m)ir)(/t + n'i^)e^('=+"'-^)^ 

77, ^ ' n' m 



^In two or three dimensions, K is replaced by the vector K which is called the reciprocal lattice 
vector in crystallography and solid state physics. 
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Prom the coefficients of all harmonics, an infinite set of homogeneous algebraic equations 
for K: 

,2 



K - X! ^mVn-m{k + (n - m)K){k + nX) = (2.25) 



In general these can be solved numerically. 

As a reference, the solution for the limiting case of all = (uniform medium) is 



where Vq is a constant. The eigenvalue condition is simply 

±k 



(2.26) 



C 



(2.27) 



which are two straight lines in the upper half oi k — uj plane. See Figure ??. 
We know, however, that the general Bragg resonance condition is 

k + nK^-k, i.e., k ^ ±K/2,±K,±3K/2, ... 



(2.28) 



Let us consider the Bragg resonance near k — NK — —k or k — NK/2, for small 
Cm = 0(e). 

We first take n = —N, 

,2 



^-{k- NKf V^^ - e_Nk{k - NK)Vo 
- J2 emV-N-m{k + i-N - m)K){k -NK)^0 



(2.29) 



and then take n — 0, 



UJ 



Vo - eNV-N{k - NK)k 



emV-rn{k - mK)k = 



(2.30) 



For any n^N, (2.25) tells us that K = 0(6)^0- Hence in (2.29) and (2.30), we can 
ignore most terms except V-n and Vq, so that 

,2 



(k - NKy 



V-N - e-Nk(k - NK)Vq = 



(2.31) 
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and 

"2 

-k'^ Vo-eN{k-NK)kV-N^O (2.32) 



C2 

Let us now shift the origin: 
k^^ + k', ^ = ^ + ^, where A;' ~ 0(e_^) < 1, ^ ~ 0(e_^) < 1. (2.33) 

By keeping only terms of the first order k', cu' we get 

N'^K'^ / uj' \ 

e^N—^Vo + NK { - + k')V.N = 



4 " ' \C 



NK{—-k']VQ^ eN——y-N = 



C 

Vanishing of the coefficient determinant gives the eigenvalue condition 

'^)' -k"-Y^\eN?N^K^ (2.34) 

which is represented by two branches of hyperbola in the k' vs. uj' plane. Use is made 
of the fact that e_jv = e^. Since 



If, 

^ > \\eN\NK (2.36) 
k' is real; the wave is propagating. If 

^ < \\eM\NK (2.37) 

k' is imaginary; the wave is evanescent. The hyperbolic branches appear as discontinu- 
ities near the point of Bragg resonance as shown in Figure 2.2. This narrow band of 
frequencies in the cu/C vs. k plane is called the forbidden band in solid state physics. 
We can assign indices to the pieces of the dispersion curve : u!i{k) for the branch in 

< < K/2, uj2{k) for K/2<k< K, uj^{k) for X < < 3X/2 etc. 
Remark: If we replace khy k + MK, (2.2) becomes 

U{x) = e^^^+^^^^ K(a;)e^"^^ = e^^" J] K(^)e^^"+^^^^ 








Figure by MIT OCW. 

Figure 2: Dispersion relation and bandgaps: Top: Correction near the first Bragg 
resonance. Bottom: Dispersion curves. 
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or 

n' 

We can rename Vn'-M — Wn', without changing the solution. Also(2.25) becomes 

Vn-J2 ^rnVn-m{k + MK + {n- m)K){k + M K + uK) ^ 

m 

or 

Vn>-M - ^ e^Vn'-M-m{k + {u' - m)K){k + n'K) = 

m 

or 

Wn' - J2 emWn'-mik + {u' - m)K){k + n'K) = 

m 

Thus the eigen- value problem remains the same. 

It is therefore possible to shift all the branches a;„(/c) vs k and plot them in —K/2 < 
k < K/2, which is called the first Brillouin zone. 



— -{k + MK + nKf 



}r,-{k + n'Kf 



yr,-{k + n'Kf 



3 Scattering by a finite layer of periodic laminate 

Let us use a one-dimensional example to describe the phenomenon. First, since many 
scatters must be involved in order for this phenomenon to be appreciable, the total region 
of disturbances must be much greater than the typical wavelength. The perturbation 
method of multiple scales can be used. Second, since reflection is strong, incident and 
reflected wave must be allowed to be comparable in order. 

Let us consider the one-dimensional scattering of elastic waves in a rod with a slightly 
periodic elasticity, 

p = constant, E = Eo{l + eD cos Kx), (3.1) 



where D is of order unity, i.e.. 



ox 



, du 

(1 + eD cos Kx) — 
ox 



(3.2) 
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We now assume that the spatial period of inhomogeneity £ = 2t^/K and the elastic 
wavelength = 27r//c = 27^^/ Eg/ p/uj are comparable. As a consequence, wave reflection 
can be significant. 

Let us first try a naive expansion, u = Uq + eui + • • • . The crudest solution is easily 
found to be 

Uo 



A 



e''=^-^'"* + c.c.. 



where c.c. signifies the complex conjugate of the preceding term, and 



27r 

T 




Eo 271 ^ 

-, or 7 = C 




(3.3) 



(3.4) 



At the next order the governing equation is 



d f ^ dui 
dx \ ° dx 



P- 



-EoD d 
~2 dx 



<iKx 



EoD d 
~2^~dx 



AKx 



+ e 



-iKx 



dun 



dx 



_|_ Q-iKx-^ I l^^^ikx-iujt 



ikA 



* ^—ikx+iLut \ 



(3.5) 



2 " 2 " J 

Clearly, when 

K^2k + 5, 5^k, (3.6) 

some of the forcing terms on the right will be close to a natural mode cxp(±i(A;,x + ut)). 
Resonance of the reflected waves must be expected. It sufices to illustrate the response 
to one of these terms, 

Eo^-p^^Ae^^^e^'^, with (t>o^kx + ut. 

Combining homogeneous and inhomogeneous solutions and requiring that ui{Q,t) — 0, 
we find 

""^^ E,{{k + 5f-k^y 

Clearly if 5 = 0(e), eui ~ 0{e/5) and is not small compared to uq except for 5x <^ 1. 
Furthermore as x increases, Ui grows as ex. Thus when Bragg condition is satisfied, the 
reflected waves are resonated and is no longer much smaller that the incident waves in 
the distance ex = 0(1). 

Let us now focus attention on the case of Bragg resonance. To render the solution 
uniformly valid for all x, we introduce fast and slow variables in space 



(3.7) 
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To allow slight detuning from exact resonance, we assume that the incident wave fre- 
quency is a; + eo;', where euj' represents the small detuning and gives rise to a very slow 
variation in time. Therefore two time variables are needed, 



t,t — et 



(3.8) 



The following multiple scale expansion is then proposed. 



u 



— Uo{x, x;t,t} + eui{x, x;t,t) + 



(3.9) 



After making the changes 



d d d d d d 
dx dx dx^ dt dt dt 



and substituting (3.9), (3.10) into (3.2), we get 



d f ^ duo 



dx dx J dt'^ 
at 0(1). Anticipating strong but finite reflection, we take the solution to be 



(3.10) 



(3.11) 



(3.12) 



where A{xi,ti) and B{xi,ti) vary slowly in space and time. At the order 0(e) we have 



d_ 
dx 



1 - = -2^0^ + 2p- 



dx 



EoD d 



2 dx 
dA^ 
dx 
dA, 



dt^ 



dxdx 



dtdt 



-E. 



+ e 

ikx—iuit 



^ ^ dx 



d R 

kJ Jb 

dB 



■ikx—iuit 



+ C.C. 



-ia;)e^'=^-^^* + c.c. + ^(-iu;)e- 
dt dt 



ikx—iujt 



+ C.C. 



EoD d 



4 dx 

The last line can be reduced to 
EoD 



I (e^*'^^ + c.c.) ^ [^e*^^-^^* + c.c. + 56"*'=^-*'"* + c.c] | (3.13) 



4 



-3ryle 



j^2Q^kx-i<^t ^ ^ ^2j^^-ikx-i^t ^ 



2 4 ^3ikx—iLut 



+ C.C. - Sk'^Be 



-3ikx—iu!t 



+ C.C.) 
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To avoid unbounded resonance of ui, i.e., to ensure the solvability of ui, we equate 
to zero the coefficients of terms e^'^^^^-'^*) and e='=*^*^^+'^*) on the right of (3.13). The 
following equations are then obtained: 

where •>/ Eo/ po — C — uj/k denotes the phase speed. These equations govern the 
macroscale variation of the envelopes of the incident and reflected waves, and can be 
combined to give the Klein-Gordon equation 

Introducing the symbol 

^ = f.^i. (3.17) 

which has the dimension of frequency, (3.16) many be written as 

With suitable initial and boundary conditions on the macro scale, one finds the slow 
variation of these wave envelopes, hence the global behaviour of wave motion. 

As an aside, let us first try a progressive-wave solution for the envelope in an infinite 
domain: 

A = Aqc'^^-^^ (3.19) 
Physically, this amounts to a detuned wave 

u = A{x, f)e''=^-*'^* = A^ ex.p[i{k + eK)x - i{uj + eQ)t\ + * , x < , (3.20) 

Clearly 

^2 _ ^2^2 ^ ^3 21) 

which represents two branches of hyperbola. Within the frequency gap \VL\ < VIq, K is 
imaginary and the envelope wave is evanescent. This result is equivalent to the band-gap 
theory before. 
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Now back to the finite laminate < x < L which is sandwiched between two semi- 
infinte sohds {x < 0,x > L) with the same uniform elasticity Eq. Let the the incident 
wave train arriving from a; ~ — oo be slightly detuned from resonance To the left and to 
the right of the laminates, the governing equations are simply 

At + CA^^O, Bt-CB^^O, x<0, and x>L. (3.22) 

We shall assume further that B — ior x > L (the rediation condition) . Over the bars 
(3.14) and (3.15), or (3.16) hold. In order that displacement and stress and horizontal 
velocity be continuous at x — 0, L, A and B must be continuous at x — 0, L. Since the 
solutions must be of the form, 

{A,B) ^ Ao{T{x),R{x))e-'^\ 0<x<L. 

T and R are governed by 

(f]2 - nl) 
nx + - ^-^r = 0, 0<x<L. 

Several cases can be distinguished according to the sign of fl"^ — Qq-. 

Subcritical detuning: < ft < Qq . 

Let 

gc= (Qg-Q')^/' (3.23) 

^ iQC cosh Q{L -x) + n sinh Q{L - x) 

iQC cosh QL + Q sinh QL ^ ' ^ 

Qsmh.Q{L — x) 
^ iQC cosh QL + n sinh QL ' ^ ''^^^ 

On the incidence side the reflection coefficient is just -R(O) and on the transmission side 
the transmission coefficient is T[L). Clearly the dependence on L and x is monotonic. 
In the limit of L ^ oo, it is easy to find that 

^W = ^""^ flW = ^e-«^ (3.26) 
Thus all waves are locahzed in the range x < 0{1/Q). 

Supercritical detuning: fl > Qq . 



then 



and 
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Let 

PC = (fi' - (]2)i/2 (3 27) 

then the transmission and reflection coifficients are: 

_ PC cos P{L-x)-msinP{L-x) 

" PC cos PL - in sin PL ^'^'^^^ 



and 



^^^^ " PCc^sPL -insJpL ■ ^^-^^^ 



The dependence on L and S is clearly oscillatory. Thus Qq is the cut-off frequency mark- 
ing the transition of the spatial variation. For subcritical detuning complete reflection 
can occur for sufficiently large L. For super-critical detuning there can be windows of 

strong transission. 

In the special case of perfect resonance, we get from (3.24) and (3.25) that 

^ ^ A cosh^^^^ ^ ^ B isinh^^^ ^ ^ 

T{x) = — = -^fj— R{x) = — = v^fr- ■ 3.30 

^ ' Ao cosh^ ^ ^ Ao cosh%^ ^ ^ 

A similar result is known for optical waves in layered media (Yariv & Yeh, 1984, p 
197), and was also found for water waves over a wavy bed. In a laboratory experiment 
for water waves, Heathershaw(1982) installed 10 sinusoidal bars of amplitude D = 5 cm 
and wavelength 100 cm on the bottom of a long wave flume. Incident waves of length 
2Tr/k — 200 cm were sent from one side of the bar patch. On the transmission side, waves 
are essentially absorbved by breaking on a gentle beach. Sizable reflection coefficients 
were measured along many stations over the bar patch. 

As shown in Figure 3, both T{x) and R{x) decrease monotonically from x = to 
a: = L, in good agreement with the experiments of Heathershaw. Thus enough small 
bars can generate strong reflection, especially in very shallow water. 

Exercise 5.1: Bragg resonance by a corrugated river bank. 

An infinitely long river has contant depth h and contant averaged width 2a. In the 
stretch < x < L, the banks are slightly sinusoidal about the mean so that 

y = ±a±BsmKx, KB = e<^l. (3.31) 
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Figure 3: Bragg scattering of surface water waves by periodic bars. Comparison of 
theory with measurements by Heathershaw. 




Figure 4: Can wavy banks serve as a breakwater? 

See Fig. 4. Let a train of monochromatic waves be incident from x ~ — oo, 

C = -e^^'^^-^*) (3.32) 

where kh,ka = 0(1). Develop a uniformly valid linearized theory to predict Bragg 
resonance. Can the corrugated boundary be used to reflect waves as a breakwater? 
Discuss your results for various parameters that can affect the function as a breakwater. 
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4 Wave localization in a random medium 



[Ref]: Mei & Pihl Localization of nonlinear dispersive waves in weakly random media, 
Proc. Roy. Soc. Lond. 2002, 458, 119-134. 

There are numerous situations where one needs to know how waves propagate through 
a medium with random impurities: light through sky with dust particles, sound through 
water with bubbles, elastic waves through a solid with cracks, fibers, cavities, hard or soft 
grains. Sea waves over a irregular topography, etc. It is known that, for one-dimensional 
propagation, multiple scattering yields a change in the wavenumber (or phase velocity) 
as well as an amplitude attenuation, if the inhomogeneities extend over a large distance. 
These changes amount to a shift of the complex propagation constant with the real part 
corresponding to the wavenumber and the imaginary part to attenuation. In particular, 
the spatial attenuation is a distinctive feature of randomness and is effective for a broad 
range of incident wave frequencies. This is in sharp contrast to periodic inhomogeneities 
which cause strong scattering only for certain frequency bands (Bragg scattering, sec 
e.g.. Chapter 1). Phillip W. Anderson (1958) was the first to show, in the context 
of solid-state physics, that a metal conductor can behave like an insulator, if the ,ir- 
crostructure has is disordered. This phenomoena, now called Anderson localization, is 
now known to be important in classical systems also. A survey of localization in many 
types of classical waves based on linearized theories can be found in Sheng (1998). 

For weak inhomogeneties, the shift of propagation constant amounts to slow spatial 
modulations with a length scale much longer than the wavelnegth by a factor inversely 
proportional to the correlation of the fiuctuations. In this section we apply the method 
of multiple scales to introduce the theory for the simplest exampole of one-dimensional 
sound. 

We begin with the Helmhotz equation for sinusoidal waves. 



Let V{x) be a random function of x with zero mean and V{x) — > 0, for x ~ —oo. An 
incident wave train 



^k'^{l + eV{x)fU = oox < oo. 



(4.1) 



dx"^ 



,ikx 



(4.2) 
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arrives from the left-infinity where there is no disorder. What will happen, on the 
average, to waves after they enter the region of disorder? 

Consider an ensemble of random media. For each realization, the wave number now 
fluctuates about the mean k by the amount order 0(e). Since {V) = 0, we expect that, 
on the average, the wave phase is affected only by the root-mean-square, wich is of the 
order O(e^). With this guess, it is natural that slow variations described by X2 — 
will be relevant. We assume that the disorder has two characteristc scales so that 

V = V{x,X2) (4.3) 

For simplicity we shall further assume that V is stationary with respect to the short 
scale 

{V{X,X2)V{XJ ,X2)) = C,v{\x - X^\,X2) (4.4) 

where (/) denotes the ensemble average of /. 
Let us try the following expansion, 

U = C/o(x, X2) + eC/i(x, X2) + e^V2{x, ^2) + • • • (4.5) 

Subsituting (4.5) into (4.1), the following perturbation equations are found, 

^ + fc^^o = 0, (4.6) 

+ k^Ux = -2A;Vt/o, (4.7) 



dx^ 



k'U2 = {2VU, + V'Uo) , (4.8) 



dx"^ dxdx2 
The solution at the leading order is 

Uo = A{x2)e'''"' where A(0) = A. (4.9) 

At the next order the inomogeneous equation is solved by Green's function G{x, x') 
defined by 

^ + eG = 5{x-x'), (4.10) 
where G is outgoing at infinities. The solution is found in Appendix A, Chapter 2 to be 



G = -^e^'^l^-^'l (4.11) 
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(4.12) 



The solution for Ui is 

/•CO 

C/i = - / dx'G{x, x') [2k'^V{x\ X2)U{x', X2)] 



-00 
f 00 



/oo 
dxV(a;', 0:2)6^'=^' e^'^l^-^'l (4.13) 
■00 

which is random with zero mean. For the O(e^) problem, we note that 

dxdx' 8x2 ' 

2eVUi^2ikA{x2)e"''' J V{x,X2)V{x' ,X2)e'''^''~'''^e-'''^''-'''^dx' , 

k^V^Uo = k^e'''''V{x,X2)V{x,X2)A{x2). 
We now take the ensemble average of (4.14), and get 

/oo 
{V{x, X2)V{x', a;2))e*^l^-^'le-^'^(^-^')(ia;' 
-00 

- k^e''''A{x2){V\x,X2)) 

For {U2) to be solvable, we set the right-hand-side to zero, 

8A ( ik "I 

^ + Alk'' j da;'C,^(|a;-a;'|,a;2)e''=l^-^'le-^''("-"')da;'- yC^^(0,a;2)j =0 

Clearly the integral above is just a known function of X2 once the correlation function 
is prescribed. Denoting 

dx'a.(k-x'|,X2)e^'=l"-"'le-^'=(^-^')dx'- -a„(0,X2). (4.14) 
■00 ^ 

If /3 = 0, X2 < and (5 — constant, X2 > 0, then the solution is simply 

A = ^(0)e-*^^^2e-^'^^2 (4.15) 

Thus, not only the phase is changed but the amplitude decays exponentially over the 
distance 0{L) where 

L=l/(3re^ (4.16) 
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In summary, due to scattering by disorder, an apparent damping is created. The distance 
L is called the localization distance. 

For simple correlation functions, the integral for /5 can be explicitly evaluated. For 
example let 

a.(k-x'|,X2) = (72(x2)e-"l^-^'l (4.17) 
so that (7^ is the RMS amphtude of the disorder. We leave it as an exercise to show that 



f 

J — c 



so that 



2 2 Q;^ + 2/^^ ika"^ 



P = Pr + 'W. = 2kV^^^^^^^-—^^^ (4.19) 



The leading order wave is 



Uq = Aq exp < ik 



1 + 



2 a^ + Ak^ 



(4.20) 

As the RMS of the disorder increases, the wwavnumber increases, hence the wave length 

decreases. A dimensionless localization distance can be defined as 

1 l + 4k^/a^ 
I'^ioc - 2e2^2 {k/a){l + 2kya^) ^ ' 

Note that the correlation length is 0{a^^). If the waves are much longer than the 
correlation length, k/a <^ 1; kLioc increases without bound and localization is weak. If 
the waves are much shorter than the correlation length fc/a S> 1; kLio^ decreases; waves 
cannot penetrate deeply into the disordered region. 

lAP (challenge) Project : Scattering of elastic waves by random distribution of 
hard grains or cavities. 



References on Optics, Physics, etc.: 
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A. Yariv & P. Yeh. Optical Waves in Crystals. Wiley 1984 

K Inoue, & K, Ohtaka, Photonic Crystals (Physics, Fabrication and Applications) , 
Springer 2004 

J-M. Lourtioz, H Bcnisty, V. Bcrgcr, J-M Gerard, D Maystre & A A.Tchelnokov, 
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